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ABSTRACT 

We present a new method to estimate the average star formation rate per unit stellar 
mass (SSFR) of a stacked population of galaxies. We combine the spectra of 600 — 1000 
galaxies with similar stellar masses and parameterise the star formation history of 
this stacked population using a set of exponentially declining functions. The strength 
of the Hydrogen Balmer absorption line series in the rest-frame wavelength range 
3750 - 4150A is used to constrain the SSFR by comparing with a library of models 
generated using the BC03 stellar population code. Our method, based on a principal 
component analysis (PCA), can be applied in a consistent way to spectra drawn from 
local galaxy surveys and from surveys at z ~ 1, and is only weakly influenced by 
attenuation due to dust. We apply our method to galaxy samples drawn from SDSS 
and DEEP2 to study mass-dependent growth of galaxies from z ~ 1 to z ~ 0. We 
find that, (1) high mass galaxies have lower SSFRs than low mass galaxies; (2) the 
average SSFR has decreased from z = ltoz = 0bya factor of ~ 3 — 4, independent 
of galaxy mass. Additionally, at z ~ 1 our average SSFRs are a factor of 2 — 2.5 lower 
than those derived from multi-wavelength photometry using similar datasets. We then 
compute the average time (in units of the Hubble time, tn(z)) needed by galaxies of a 
given mass to form their stars at their current rate. At both z = and at z = 1, this 
timescale decreases strongly with stellar mass from values close to unity for galaxies 
with masses ~ 10 10 Mq, to more than ten for galaxies more massive than 10 u Mq. Our 
results are in good agreement with models in which AGN feedback is more efficient at 
preventing gas from cooling and forming stars in high mass galaxies. 

Key words: galaxies: evolution - galaxies: star formation 



1 INTRODUCTION 

Over the last decade, there have been a large number of 
photometric and spectroscopic surveys designed to study 
the formation and evolution of galaxies. One major conclu- 
sion of these studies has been that the epoch when massive 
galaxies formed most of their stellar mass is significantly 

* Email: chenym@mail.ihep.ac.cn 



earlier than that for low mass galaxies dHeavens et alj|2004l ; 
iThomas et all 2005h This phenomenon, popularly known as 
"down-sizing", at first sight seems at odds with the predic- 
tions of the hierarchical CDM model, in which dark matter 
halos of all masses grow through merging and accretion right 
up to the present day. The only way to reconcile the obser- 
vations with the theory, is to postulate that the growth of 
the most massive galaxies is much slower than the growth 
of their surrounding halos. 
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The most natural way to achieve this is by in- 
voking feedback processes, which prevent gas from cool- 
ing, condensing and forming star s in massive halos 



Croton et aHl2006l : lDe~Lucia et al.ll2006l ; 



llBower et al.ll2006l; 

Guo fc Wh ite 2008). A variety of different feedback mech- 
anisms have been included in numerical and semi-analytic 
models, for examp l e feedback from quasars and radio AGN 



( Silk fc Reeslll998l ; iHopkins et ail 12006): IWang et all l200rj. 



. . /ang 1 

20071 ) , supernova heatin g ( Cole et all " 2000l ; iBenson et al.l 



2003l ; )Stringer et al.ll200Sl) and heating by infalling substruc- 
ture |Dekel fc Birnboiml 2008), A physical understanding of 
feedback is still lacking, however, and it is not understood 
which, if any, of the proposed mechanisms are most impor- 
tant in regulating the growth of galaxies. It is likely that each 
mechanism will come into play at a different mass scale and 
cosmological epoch. By quantifying in detail how galaxies 
of different masses grow as a function of time, we hope to 
clarify how different galaxies form their stars and how this 
is influenced by feedback processes. 

There have been many attempts to measure the star 
formation rates of galaxy populations from the present 
day out to redshifts g r eater than 5 (se e for example 



Brinchmann et al. 



2005; lNoeske et all 



tts g r eater than o (se e tor exa mple 
20041; iBauer et al. I 120051 ; iFeulner et al.l 



2007aT l . For nearby galaxies, the Ha line 



provides the most reliable estimate of SFR, since it directly 
measures the number of ionizing photons from massive stars. 
A reasonably reliable correction for dust extinction can be 
made provided one also measures the Balmer decrement 
Hq/H/3 accounting accurately for stellar absorption. Beyond 
redshifts of ~ 0.4, the Ha redshifts out of the optical part 
of the spectrum and is no longer accessible. [O 11] equivalent 
widths, ultraviolet (UV)-optical spectral energy distribution 
(SED) fitting and infrared photometry are commonly used 
to estimate the star formation rate. Each of these indica- 
tors is subject to different disadvantages. The [O 11] equiv- 
alent width is strongly affected by dust and by metallicity. 
The UV luminosities of galaxies are also strongly affected 
by dust, while the infrared only provides a direct measure of 
SFR if one has measurements across the thermal peak of the 
spectrum. This has not been the case for many of the Spitzer 
surveys aimed at quantifying star formation in high-redshift 
galaxies. All the indicators discussed above may also be con- 
taminated by AGN emission if there is an actively accreting 
black hole in the galaxy. 

In this paper, we develop an approach to estimate the 
amount of recent star formation experienced by a popula- 
tion of galaxies. Our method is based on the Balmer ab- 
sorption lines located in the rest-frame wavelength range 
3750 — 4150A of the galaxy spectrum. The advantages of 
our method are the following: 

(f) The Balmer absorption lines are weakly influenced 
by dust attenuation and AGN contamination compared with 
the indicators discussed above. 

(2) The wavelength range spanned by the Balmer ab- 
sorption lines is accessible out to redshifts greater than 1, 
even in optical spectra. This means that the method can be 
applied in a consistent manner to both low and high redshift 
(z ~ f) samples. 

(3) By stacking a large number of galaxies, we estimate 
SFRs for a complete sample of galaxies in a given stellar 
mass range, including those galaxies that are forming stars 



weakly or not at all. These galaxies are often excluded when 
SFRs are measured for individual objects. 

We apply our method to a large sample of galaxy spec- 
tra from the Sloan Digital Sky Survey (SDSS) and the 
DEEP2 redshift survey to study the evolution of galaxies 
from z = 1 to z = 0. This redshift interval accounts for 
roughly half the age of the universe. Our study thus ad- 
dresses the final stage of galaxy build-up in the Universe. 

This paper is arranged as follows. In §2, we introduce 
the SDSS and DEEP2 samples used in our studies. The 
method to estimate the amount of recent star formation in 
our galaxies is developed in §3. We apply the method to 
the SDSS and DEEP2 samples, and present our results in 
§4 and §5, respectively. A discussion of the results is given 
in §6. §7 contains the summary of the paper. We use the 
cosmological parameters Ho = 70 km s" 1 Mpc^ 1 , Qm = 0.3 
and Qa = 0.7 throughout this paper. 



2 SAMPLE SELECTION 
2.1 The SDSS galaxy sample 

The lo w redshift galaxy sample comes from data release 4 
(DR4, lAdelman-McCarthv et al.|[200r3 ) of the SDSS, which 
contains more than 550000 spectra in the primary redshift 
range ^ 2 ^ 0.3. Stellar masses of the galaxies are 
obtained by fittin g a suite of BC03 models (described in 
ISalim et ail 120071 ) to the SDSS m odel magnitudes. These 
m asses are not id e ntical to those of Kauffm ann et alj (|2003l ) 
or iGallazzi et ail (|2005l ). who use spectral indices, such as 
the D4000 and H<5a to constrain the mass-to-light ratios of 
the galaxies. However, the differences between the mass es- 
timates are very small, the median offset of log M» is only 
0.01. In this work we choose to use the photometric masses in 
order to maintain consistency between our high and low red- 
shift galaxy samples. The masses assume the universal initial 
mass function (IMF) as parametrized by iKroupal l|200ll ). 

The aim of the first step of our work is to compare the 
SSFRs derived from the Balmer series with the SSFRs of 
Brinchmann et al. (2004, hereafter B04) derived from neb- 
ular emission lines. We begin with a subset of galaxies with 
14.5 < r < 17.77 and 0.005 < z < 0.22 drawn from the 
DR4 spectroscopic sample. We divide galaxies with stellar 
masses in the range of 10 9 ~ 1O 12 M0 into 6 mass bins. The 
mass interval is set to A log M, — 0.5. The SDSS spectra are 
obtained through a 3" circular fibre aperture, and therefore 
sample primarily the inner regions of galaxies. To minimize 
the effects of this so-called "aperture bias" , we select galax- 
ies with 0.9 < z/zmax < 1, where z ma x is the highest redshift 
at which the galaxy in question would pass the sample selec- 
tion criteria. In Table 1, we list the mass range and median 
redshift of galaxies in each mass bin. 



2.2 The high redshift galaxy sample 

In order to obtain accurate galaxy stellar masses at redshifts 
approaching unity, photometry in both optical and near- 
IR passbands is required. Our high redshift galaxy sample 
is selected from fields covered by both the DEEP2 galaxy 
redshift survey and the Palomar Observatory Wide Infrared 
Survey (POWIR). 
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Table 1. The mass ranges and median redshifts of galaxies in the 6 mass bins of the SDSS sample and 4 mass bins of the DEEP2 sample. 



log M*/M = 


9.0-9.5 9.5-10.0 10.0-10.5 


10.5 - 11.0 


11.0 - 11.5 


11.5 - 12.0 


z(SDSS) 


0.05 0.07 0.09 


0.14 


0.17 


0.20 


z(DEEP2) 


0.85 


0.85 


0.85 


0.85 



2.2.1 The DEEP2 Survey 

The DEEP2 Galaxy Redshift Survey (|Davis et al.l 120031. 
Faber et al. in p rep) utilizes the DEIMOS spectrograph 
l|Faber et al.ll2003l ) on the KECK II telescope. Targets for 
the spectroscopic sample were selected from BRI photome- 
try taken with the 12k x 8k mosaic camera on the Canada- 
France-Hawaii Telescope (CFHT). The images have a lim- 
iting magnitude of 7?ab ~ 25.5. As the _R-band images 
have the highest signal-to-noise ratio (S/N) of all the CFHT 
bands, they were used to select the DEEP2 targets. The 
CFHT imaging covers four widely-separated regions, with a 
total area of 3.5 deg 2 . In fields 2-4, the spectroscopic sample 
is preselected using ( B — R) and (R — I) colors to eliminate 
objects with z < 0.7 |Davis et al.l[20 03) . Color and apparent 
magnitude cuts were also applied to objects in the first field, 
the Extended Groth Strip (ECS) , but these were designed to 
downwei ght low redshift galax ies rather than eliminate them 
entirely (|Willmer et al .1 120061 ). The third data release of the 
DEEP2 survejQ contains spectra of about 50000 galaxies in 
the magnitude range 18.5 ^ Rab ^ 24.1. The spectra have 
a resolution of about R ~ 5000. 



2.2.2 The Palomar Observatory Infrared Survey 

7f s -band photometry was obtained in portions of all four 
fields targeted by the DEEP2 galaxy redshift survey us- 
ing the Wide Field Infrared Camera (WIRC, Wilson et al. 
2003 ) on the 5m Hale T elescope at the Palomar Observa- 
tory. iBundv et alj |2006l ) mapped the central third of fields 
2-4 using contiguously spaced pointings tiled in a 3 x 5 pat- 
tern. About 85% of the pointings have a depth greater than 
Kab = 22 and the imaging covers a total of 0.9 deg 2 . In field 
1 (the Extended Groth Strip), the iG-band imaging covers 
0.7 deg 2 to a typical depth of greater than Ka b = 22.5. 

W e use the stellar masses calculated by IBundv et al.l 
(2006). They estimate the K s -band mass-to- light ratios 
for the galaxies in this sample by comparing the opti- 
cal to infrared spectral energy distributions (SEDs) to a 
grid of 13400 model SED s generated using the BC03 code 
i|Bruzual fc Chariot! 12051 ) . The typical error in the stellar 
mass estimates is about 0.2dex. A IChabrierl (|2003l ) IMF 
is used: the systematic difference between the SDSS and 
DEEP2 stellar masses caused by the different choice of IMF 
is about 0.05dex. We correct for this by adding 0.05dex to 
the logarithm of the DEEP2 galaxy masses. 



1 http://deep.berkeley.edu/DR3/ 
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2.2.3 Our high redshift galaxy sample 

The high redshift galaxy sample used in this work combines 
data from the DEEP2 and POWIR survey. The criteria used 
to select the galaxies are the following: 

(i) The redshift range is 0.75 ^ z ^ 1. This ensures that 
all the galaxy spectra fully cover the rest-frame wavelength 
range required to measure the spectral indices analyzed in 
this paper (3750 - 4150A) (Note that the DEEP2 spectra 
span the wavelength range 6500 - 9100A). 

(ii) The i?-band magnitude range is 18.5 ^ Rab ^ 24.1. 
Rab = 24.1 is the limit of the DEEP2 spectroscopic sample. 
By co mparing with a fain ter sample with photometric red- 
shifts, IBundv et al.l l|2006l ) concluded that the spectroscopic 
sample is essentially complete down to M* ^ 1O 1O M0. As 
described in §4.1 we will correct for any residual incomplete- 
ness by applying a 1/V ma x weight to each galaxy. 

(iii) The JG-band magnitude limit is Kab ^ 22. The 
POWIR _K\,-band survey covers different areas to differ- 
ent depths. Kab ^ 22 is chosen to maximize the sky 
area covered, but also maintain a high level of complete- 
ness for galaxies with stellar masses logA^/M© ^ 10.0. 
This i^s-band limit eliminates around 15% sources with 
10 < log M»/M e s: 10.2 from the Rab < 24.1 sample. We 
show in §4.1 that this does not affect our results. 

In summary, our DEEP2 sample has a redshift range 
z ~ 0.75 - 1.0, with 18.5 < Rab < 24.1 and Kab «S 22. It 
covers a total area of 1.6 deg 2 and the total source number is 
about 3000. The galaxies have a mean and median redshift 
of 0.87 and 0.85 respectively. 



3 THE SPECTRAL INDICES 

The 4000A wavelength region of the galaxy spectrum con- 
tains abundant information for constraining the recent star 
formation histories of galaxies. Two traditionally-used in- 
dices, D4000 and H5a, are located in this region. Impor- 
tantly, this rest-frame wavelength region is usually included 
in both low and intermediate redshift galaxy surveys. 

PCA is a standard multivariate analysis technique, de- 
signed to identify correlations in large datasets. Using PCA, 
Wild et al. (2007, hereafter W07) developed a set of new 
high signal-to-noise ratio spectral indicators located in the 
rest-frame wavelength range 3750 — 4150A. At wavelengths 
around 4000A galaxy spectra vary in both spectral shape 
and strength of the Hydrogen Balmer absorption lines. 
These two are inversely correlated: younger galaxies have 
stronger Balmer absorption lines and weaker 4000A break 
strengths. The third axis of variation is the CaII(H&K) 
lines, which are related to both the age and metallicity of 
the galaxy. 

Because the DEEP2 spectra are not flux calibrated, the 
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continua contain no useable information. In this work, we 
do not use the spectral indicators given by W07 but closely 
follow their method to create a similar set of indicators de- 
signed to work on high pass-filtered spectra. Because the 
continuum is removed by the filter process, we expect our 
eigenspectra to reveal the Balmer absorption lines as a pri- 
mary axis of variation and CaII(H&K) lines as the secondary 
axis. We refer the reader to W07 for more details on the PCA 
method. Here we present the dataset and the methods used 
to obtain the indices and the resulting eigenspectra. 



3.1 Input model data set for PCA 

Our input data set for the creation of the PCA eigenspectra 
is a set of model spectra g enerated using the BC03 stellar 
population synthesis code l|Bruzual fc C hariot 2003) . The 
model libra ry is similar to th at used in Kauffmann et all 
(|2003h and ISalim et all (120051 ). although with a more re- 
stricted parameter range. 6641 model galaxies are selected 
at random from the parent library according to the following 
criteria: 



(i) The time t{ olnl when the galaxy begins to form its stars 
is distributed uniformly between and 5.7 Gyr after the Big 
Bang (the age of the universe is assumed to be 13.7 Gyr). 

(ii) The model galaxies have exponentially declining star 
formation histories SFR oc exp(— t/r) with r distributed 
uniformally between 1 ^ r 1.4 Gyr. 

(iii) Top-hat bursts are superimposed on these continuous 
models. Two parameters describe the bursts: f burst, the frac- 
tion of the total stellar mass formed in bursts, is distributed 
logarithmically between 0.0 and 0.1; tt, U rst, the duration of 
the burst, is distributed uniformly between 0.03 and 0.3 Gyr. 
During the burst, stars form at constant rate. Bursts occur 
with equal probability at all times after if orm and the prob- 
ability is set so that 50% of the galaxies in the library have 
experienced a burst over the past 2 Gyr. 

(iv) The metallicity is distributed linearly in the range 
0.5 ^ Z ^ 2Zq; no metallicity evolution is included. 

These specific criteria (e.g. the distribution of tf orm ) 
are chosen so that the model galaxies span a reasonably 
wide and even range in mean stellar age. In this case, the 
Balmer absorption lines and CaII(H+K) absorption lines 
will be isolated in the output eigenspectra, and the Balmer 
absorption lines will dominate the first eigenspectrum and 
CaII(H+K) absorption lines the second eigenspectrum. We 
note that although changes to the input library will alter the 
resulting eigenspectra, our calibration using BC03 models in 
Section 5 ensures that we will recover the same SSFRs. 

We convolve the model spectra to have a velocity dis- 
persion equal to 150 km/s; this is the median value of the ve- 
locity dispersion of the galaxies in our samples. Each model 
spectrum is then normalised by the mean flux in the rest- 
frame wavelength range 4000 — 4080 A, where the spectrum 
is free from strong absorption lines. 

We smooth the spectrum F\ with a Gaussian kernel 
W(X). This yields the low-pass spectral component 



Fx 



where 



/ d\'W(\ - X')Ivar(X')F x , 
J d\'W(X- X')Ivar(X') 



(1) 






TV 








3700 3800 3900 4000 4100 4200 

wavelength (A) 

Figure 1. Examples of the filter process applied to the model 
spectra. Black: the original model spectrum. Blue: the low-pass 
filter component. Red: the resulting continuum filtered spectrum, 
shifted upwards for easy comparison with the original spectrum. 
Green: the zero point. Models have been selected to cover the full 
range of galaxy ages: galaxy age increases from top to bottom. 



W(X — A') 



27TCT 



exp 



A — A' 



(2) 



and a = 32A. Ivar(X') is the inverse variance matrix of 
the spectrum. For the model spectrum, it is equal to 1 at 
all wavelength points. In the filtering process the strongest 
absorption lines are masked by setting Ivar(X') = 0. We 
mask 20A centred on A4103, 3889, 3835; 16A centred on 
A3798; 12A centred on A3770, 3750. The sizes of these masks 
were chosen to ensure that the absorption lines were not 
diluted in the filter process. Small changes to the mask sizes 
make no difference to the final results. 

The high-pass component is then obtained as 



F 



Fx - Ft 



(3) 



We choose to subtract the low-pass component, rather than 
divide by it as is commonly done, because we wish to retain 
the luminosity weighting of the real galaxies when they are 
stacked (see the following Section). Figure 1 shows three 
examples of the filtering process. 



3.2 Creating the 'eigenspectra' 

The model library is dominated by old galaxies. To ensure 
the dominance of the Balmer absorption lines in the first 
eigenspectrum, we weight each model spectrum by the re- 
ciprocal value of its 4000 A break strength to reduce the 
impact of galaxies without strong Balmer absorption lines. 
The mean spectrum of the input data set is then calculated 
and subtracted, and the PCA is run to produce the eigen- 
spectra. 

Figure 2 presents the mean spectrum and first two 
eigenspectra of our input model galaxies. As required, 
the first eigenspectrum shows the Balmer series. The 
CaII(H&K) absorption lines are isolated in the second eigen- 
spectrum. In the following section the first 10 eigenspectra 
are used to fit the SDSS and DEEP2 spectra. The coef- 
ficient of the first eigenspectrum (PCI, see eq.14), which 
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Figure 2. Top: the mean spectrum of the continuum filtered in- 
put data set. Middle: the first eigenspectrum in which the entire 
Balmer absorption line series is visible. Bottom: the second eigen- 
spectrum, which shows the CaII(H&K) absorption lines. 



represents the strength of Balmer absorption lines, is used 
to derive SSFRs. 



4 DATA ANALYSIS 

In this section, we describe how we apply our method 
to real data. O u r me thodology is similar to that of 
iGlazebrook et al l (|2003h . who stacked a large number of 
SDSS spectra to create a "cosmic average" spectrum of all 
the stars in the Universe. The main difference is that we 
carry out the cosmic averaging in bins of stellar mass. 
The steps of the analysis are the following: 

• The galaxy spectra are corrected for foreground Galac- 
tic extinction and shifted to the restframe. 

• The DEEP2 spectral resolution is about 3 times higher 
than that of the SDSS and the BC03 population synthesis 
models. In order to analyse SDSS and DEEP2 spectra in an 
equivalent manner, the DEEP2 spectra are convolved with a 
Gaussian kernel to match the resolution of the SDSS spectra. 

• The spectra are rebinned to match the model wave- 
length binning. 

• The spectra are converted into units of luminosity (see 
below for details). 

• The continuum of the spectra is subtracted using the 
same filtering method used for the model spectra. An extra 
mask including 14A centered on Ne IIIA3783 line is added. 

All the filtered spectra in the same mass bin are stacked 
using the following weighting scheme: 



■p(A) = 



5Z"=i W(i) x Ivar(i, A) x L(i,X) 
J27=i w ( i ) x Ivar(i,X) 



(4) 



where the sum is over i = 1,2, ...,n galaxies with individ- 
ual luminosities L(i,X) in units of ergs _1 A . Ivar(i,X) is 
described below, and accounts for bad pixels and night sky 
linefl W(i) provides a weight for each source to correct for 
missing galaxies due to survey incompleteness. 



For the DEEP2 spectra we make use of the code 



For the flux calibrated SDSS spectra we convert to lumi- 
nosity by scaling by the luminosity distance of the galaxies in 
the standard way. We set Ivar(i, A) = for bad pixels (iden- 
tified in the SDSS mask array), otherwise Ivar(i,X) — 1. 
The redshift intervals have been selected so that the sam- 
ples are complete in mass and we therefore set W(i) = 1 for 
all galaxies. 

For the DEEP2 galaxies, the spectra are not flux cal- 
ibrated and cannot be converted directly into luminosity. 
In order to get L(i,X), we first normalise each individual 
spectrum over all wavelengths to its mean flux between 
4000 — 4080A. This normalised spectrum is then scaled ac- 
cording to the optical photometry using the known spec- 
troscopic redshift and optical colour information. The K- 
correction is derived from the best -fitting SEP template out - 
put by the KCORRECT code (|Blanton fc Roweisl I2007T I . 
The inverse variance matrix Ivar includes only variations 
due to sky intensity, not object intensity 2 . This ensures that 
bad pixels are downweighted while equal weight is given to 
each spectrum. 

Setting Ivar(i, A) = 1 and down-weighting the noisy 
pixels give similar results for SDSS. But for DEEP2, these 
two methods give very different results. In order to get the 
same S/N of the composite spectra, we need to stack many 
more galaxies if we set Ivar(i,X) = 1 than if we down- 
weight the pixels effected by sky lines in the DEEP2 spec- 
tra. This is because local galaxy spectra (SDSS) have much 
much higher S/N than high redshift galaxies (DEEP2) on 
average. The reason that individual pixel noise arrays are 
not used in the Ivar array, is that this would overall down- 
weight fainter galaxies, possibly introducing biases that are 
difficult to quantify into the composite spectra. 



4.1 DEEP2 completeness correction 

An issue with every data set is the selection of weights to 
correct for missing galaxies. The weights we use in this study 
take into account three factors: 

• The target selection strategies of DEEP2. 

• The redshift success rate. 

• The magnitude limits. 

The weight W(i) for DEEP2 galaxy i is expressed as: 



W(i) 



\'' 
i in 



(5) 



where m accounts for incompleteness resulting from the 
DEEP2 color selection and redshift success rate. V^a* ac- 
counts for the fact that faint galaxies are not detected 
throughout the entire survey volume in a magnitude-limited 
survey. 

To calculate Kj we closely follow the method described 
in lWillmer et al.l l|2006l ). but add an extra dimension to the 
reference data cube, which is the iCVband magnitude. For 
each galaxy in the photometric catalogue, we select all the 
galaxies from the redshift catalogue sharing the same bin 
in the (B — R)/(R — I) / Rab/ K s space, count the number 
of objects with failed redshifts (Nf); the number of galaxies 



COADD_SPECTRA.PRO available from the DEEP2 Team to identify 
bad pixels and separate sky from continuum noise. 
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with z < zi (N ZI ); the number of galaxies with z > z u (N Zu ); 
and the number of galaxies with successful redshifts within 
the "legal" redshift range z t -z u (N z ). For field 1 (the EGS), 
zi=0.2, Zu=lA; for fields 2-4, z t =0.7, z u = lA. The proba- 
bility that each galaxy in the photometric catalogue has a 
redshift in the legal range is estimated. For galaxies with 
high-quality redshifts, the probability P(zi ^ z ^ z u ) — l 
when zi $J z z u and P(zi ^ 2 ^ z u )—0 if z < zi or z > z u . 
The estimation of the probability for unobserved sources is 
based on the so-called "optimal" weighting model (see §2 of 
Willmer et al. 2006), which assumes that the failed redshifts 
of red galaxies have the same distribution as the successful 
ones, while blue galaxies with failed redshifts lie beyond the 
redshift limit z u . Namely, for the red galaxies: 

p{zi < * < *•> = n. + n: 1+n ^ (e) 

for the blue galaxies 



MF: z=0.75-1.0 



P(z t < z < z u ) = 



N, 



+ N ZI + N Zu + N f 

Finally, m is calculated by summing over all galaxies j in 
the photometric catalogue with the same (B — R)/(R — 
I)/Rab/K s value as galaxy i, the probabilities that the red- 
shifts of galaxies are within the legal limits 



N z 



(8) 



In the case of EGS, a final correction is applied to m to 
account for the different sampling strategy, which includes 
low-redshift (z < 0.7) galaxies but downweights them so 
that they do not dominate the sample. The correction fac- 
tor (/ m ) depends on the location of the galaxy in (B — R) 
versus (R — 7) space and its apparent magnitude. / m has 
been computed for aroun d a quarter of the EGS galaxies 
from lWillmer et all (|2006l ) for the purpose of measuring the 
galaxy luminosity function. We estimated f % m for the remain- 
ing galaxies in the EGS data to be the average f m of all the 
galaxies within the same bin in the (B — R)/(R — I)/Rab- 
The probability that a galaxy will be placed on an EGS 
mask is given by 



P(mask) = 0.33 + 0.43P ga i/„ 



(9) 



where P ga i, given in the DEEP2 catalogue, is the probability 
that an object is a galaxy based on its magnitude, color and 
size. For EGS galaxies, 

** ~ iV z P(mask) (W) 

where j includes all the galaxies in the photometric cata- 
logue which are within (B — R)/{R — I) I RabI Ks spa ce. 

The calculation of V^ ax follows ISchmidt] l|l968l ). and 
provides a simple way to account for the P-band limit of 
the sample. We define 



dV 

dz 



dz 



(11) 



where df2 is the solid angle covered by the sample, 4*- is 
the comoving volume element and z m i n ,i and z ma x,i are the 
low and high redshift limits within which galaxy i can be 
detected given the P-band limit of the survey. They are given 
by: 




12.0 



Figure 3. The stellar mass function for galaxies in the redshift 
range 0.75 — 1. The black line is our result, the pink line is from 
Bundy et al. (2006). 



Zmax.i = mm{z' max ,z(Mi,m u )}, 
Zmin.i = max{z| ]lill ,z(M i ,m ! )}, 



(12) 
(13) 



where z' min and z' max are the redshift limits of the sample, 
z(Mi,m u ) is the redshift above which the galaxy would 
be fainter than the P-band magnitude limit of 24.1 and 
z(Mi,mi) is the redshift below which the galaxy becomes 
brighter than Pab = mi = 18.5. We use the best-fitting SED 
templates output by KCORRECT to calculate z(Mi,m u ) 
and z(Mi,mi). Because we are analyzing galaxies in a rel- 
atively narrow redshift slice, no evolutionary correction is 
applied. 

In order to check that we recover the correct galaxy 
weights, we calculate the stellar mass function for the 
redshift bin 0.75 — 1 and compare it to the result from 
I Bundv et alj {2006). Figure 3 shows that the two stellar 
mass functions are almost identical. 

In §2.2.3, we noted that some sources are missing in 
our log M*/M Q = 10 - 10.5 DEEP2 mass bin due to the 
additional Pab = 22 magnitude limit. To test whether this 
incompleteness affects our final results, we repeat the entire 
analysis for a limiting magnitude Pab = 22.5. This limit 
corresponds to a mass completeness limit of about 10 10 Mq. 
We find that our measured Balmer absorption line strengths 
in the lowest mass bin of DEEP2 remain the same to within 
the errors. If we were to increase the P s -band magnitude 
limit to this fainter value for the whole sample, the area 
covered by the sample would decrease and this would affect 
our sample statistics, particularly for our highest mass bins. 



4.2 Fitting composite spectra 

The PCA pr ojection algorithm emplo yed here accounts for 
pixel errors (|Connollv fc Szalavl Il999h and allows the nor- 
malisation of the spectrum as a free parameter (G. Lemson, 
private communication, see W07). The fitting result is a lin- 
ear combination of the eigenspectra, 



ap (A) =N x 



(14) 
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Table 2. Measured values of iib a i mcr and the SSFRs derived from them for each mass bin of SDSS and DEEP2. 'S' represents SDSS, 
'D' represents DEEP2. 
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Figure 4. The composite spectra (black) and PCA fitted models (red, by eq.14) for all mass bins. The overall normalisation (TV) has 
been divided out to make comparison easier. 



where JV is the normalization, M represents the mean spec- 
trum (Figure 2, top). Pd is the amplitude of the i— th eigen- 
spectrum Ei, representing the amount of the i— th eigenspec- 
trum present in the residual spectrum. These amplitudes are 
the new spectral indices which are calculated from the pro- 
jection of each composite spectrum onto the eigenspectra. 
We note that principal component amplitudes can be nega- 
tive, for example, negative PCI means the strength of the 
Balmer absorption lines in the composite spectrum is weaker 
than that in the mean spectra. 751 represents the first eigen- 
spectrum (Figure 2, middle), etc. In the current work, we are 
interested in the first index, which measures the strength of 
the Balmer series. We refer to it as 'i/baimer' from now on. 
We would like to stress that because we stack all spectra in 
each mass bin, the value of ffbaimor relates to the "cosmic 



average" spectrum, which includes old galaxies with weak 
Balmer absorption lines. 

Figure 4 shows the normalized composite spectra 
(£comp(A)/iV ', black) and the the best-fit results (red). Both 
the strength of the Balmer absorption and the nebular emis- 
sion lines in these spectra reveal that for the same mass bin, 
DEEP2 galaxies have much younger populations than SDSS 
on average. The values of -Hbaimer are listed in Table 2 for 
each mass bin of SDSS and DEEP2. Although statistical er- 
rors on the principal component amplitudes are calculated 
during the projection of the eigenspectra, they are insignif- 
icant with respect to the errors caused by the sample size. 
We estimate the sample variance errors using a bootstrap 
technique, by sampling with replacement the same number 
of galaxies as are in the stack, and recalculating i^baimer 100 



© 0000 RAS, MNRAS 000, 000-000 



8 Chen et al. 



times. The errors are given by the root-mean-square differ- 
ence between the -Hbaimer values calculated for each boot- 
strap sample and the mean value. 



5 ESTIMATING SPECIFIC SFR 

When stacking many hundreds of galaxy spectra, it is rea- 
sonable to assume that stochastic fluctuations in the star for- 
mation histories of individual galaxies will average out, and 
that the star formation history of the composite "galaxy" 
can be approximated by a smooth function. In this section, 
we work with a grid of exponentially declining star forma- 
tion models constructed from the BC03 population synthesis 
code. We create 10 6 different models by varying the follow- 
ing parameters: 

(1) The formation time of a galaxy (tf orm ) is defined as 
the time since the Big Bang when the galaxy first begins to 
form its stars and is uniformly distributed between and 10 
Gyr. 

(2) The star formation rate is modelled as SFR oc 
exp(— t/r), with an e-folding timescale r uniformly dis- 
tributed in the range 0.01 — 10 Gyr. 

(3) We calculate three sets of models with metallicities 
Z = 0.02 = Z Q , Z = 0.00 8 and Z = 0.004. According to the 
mass -metallicity relation (jGallazzi et alj|2005l ; IPanter et al.l 
2008), the average global metallicity is near solar for the 
galaxies with stellar masses greater than 1O 1O M0. For galax- 
ies with stellar masses smaller than 10 10 Mq, most galax- 
ies have sub-solar metallicity, but the metallicity increases 
rapidly with stellar mass. In the following, we will assume 
the observed mass-metallicity relation and use solar metal- 
licity models for mass bins with M* 1O 1O M0, models with 
Z = 0.008 for the mass bin with 9.5 < log M»/M < 10, and 
Z = 0.004 models for the mass bin with 9 < log M*/M < 
9.5. 

(4) We consider models with and without dust. The 
effects of dust attenuation on the spectral properties are 
computed acco r ding t o the simple two-component model of 
ICharlot fc Fail (|2000h . The total effective V-band optical 
depth tv of the models is set to be a Gaussian distribution 
with peak at 1.6 and standard deviation of about 0.8. If 
tv < 0, we set it to zero. This distribution is similar to the 
SDSS DR4 star forming galaxies with stellar masses in the 
range of 10 9 — 1O 12 M0. The parameter /x, the fraction of the 
dust optical depth contributed by the 'ambient' interstellar 
medium, is set to be 0.3. We will show later that our results 
are not greatly affected by dust. 

We extract the spectrum for each of our model galax- 
ies at an age a<?e m odci = tu{z) — tform, where £h(z) is the 
age of universe at redshift z, and z corresponds to the me- 
dian redshifts of our SDSS and DEEP2 samples. They are 
listed in Table 1. After broadening the spectra to a veloc- 
ity dispersion of 150km/s and filtering out the continua, we 
measure the values of -Hbaimer for each of the model spec- 
tra. The parameters of the models that best fit the -Hbaimer 
values measured for the composite spectra, are obtained by 
weighting each model in the library by the probability func- 
tion exp(— X 2 /2), and then binning the probabilities as a 
function of the given parameter value (see Appendix A of 
Kauffmann et al. 2003 for a detailed description). The most 
probable value of a model parameter for a particular com- 




2 4 6 8 10 

T (Gyr) 



Figure 5. Probability contours of tf orm versus r for galaxies 
in different stellar mass bins. They represent the 68% confi- 
dence regions. Results are shown for SDSS (light grey areas) 
and for and DEEP2 (grey areas). The labels 'S' and 'D' in this 
plot represent SDSS and DEEP2, respectively. 'SI' means the 
first mass bin of SDSS, with log M»/M Q = 9.0 - 9.5; 'S2' with 
log M*/Mq = 9.5 - 10.0; 'S3' with log M,/Mq = 10.0 - 10.5; 
'S4' with log M»/M Q = 10.5 - 11.0; 'S5' with log M,/Mq = 
11.0 - 11.5; 'S6' with logAf„/M = 11.5 - 12.0; 'DP with 
log M„/Mq = 10.0 - 10.5; 'D2' with log M*/Mq = 10.5 - 11.0; 
'D3' with logA/./Afo = 11.0 - 11.5; 'D4' with LogM*/Af Q = 
11.5 - 12.0. 

posite spectrum can be taken as the peak of this distribu- 
tion; the most typical value is its median. In the following 
sections, we use the median of the probability distribution 
as our adopted best estimate. 

Figure 5 shows our best-fit results in the space of r ver- 
sus tform- The shaded areas denote 68% confidence regions 
with | .Hbaimer (m) — -Hbaimer (d) | < If, where Hbaimer (m) is 
the value of -Hbaimer measured from the model galaxy spec- 
trum. -Hbaimer (d) is measured from the composite spectrum, 
a is its error. Light grey regions correspond to the 6 stel- 
lar mass bins for the SDSS sample, while dark grey regions 
correspond to the 4 DEEP2 mass bins. The central black 
lines show the best-fit results. As shown in Figure 5, there 
is a strong degeneracy between r and tform i so the l -0 " con- 
tours appear as strips, rather than compact regions in the 
plot. However, it is clear that the mass bins separate clearly 
in T-if orm space, indicating different mass galaxies have had 
significantly different recent star formation histories. 

Rather than using our model grid to estimate the pa- 
rameters tform and r, which we have shown are degener- 
ate, we now turn to a parameter that directly character- 
izes the recent star formation histories of the galaxies that 
make up our composite spectra. We remind the reader that 
the Balmer absorption lines are most sensitive to stars that 
have formed over the past 10 8 — 10 9 years. In the spectra 
of individual galaxies, the Balmer absorption lines have of- 
ten been used as a post-starburst indicator, because they 
remain strong for nearly a Gyr after the burst has already 
ended. For our composite spectra, however, bursts are not 
of relevance and our approach in this paper will be to use 
fbaimer as a measure of the average specific star formation 
rate (SSFR= (SFR)/(M„)) of the galaxies in the stack. 

We use the Balmer absorption lines as our SSFR indi- 
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Figure 6. The probability distributions of the SSFRs for SDSS 
galaxies. Lines with different colors represent different mass bins. 
The SSFRs decrease with increasing mass. Solid lines show SSFRs 
estimated using BC03 models without dust, while dashcd-lincs 
show the results using models with dust. 



DEEP2 SSFR 



0.0 















Log M./Mo 


no dust I 




D1: 


10-10.5 


dust 




D2: 


10.5-11 


D3 : 












D3: 


11 -1 1.5 






D4: 


1 1.5-12 










/ ; \ \ ! \ \ 01 : 






04 













-12 



-11 -10 

Log (SSFR /yr"') 

Figure 7. Same as Figure 6, but for DEEP2. 



cator. These are mainly contributed by A stars (see Figure 
1 of W07) , which have a lifetime of several hundred million 
years. So long as the star formation rates of our compos- 
ite galaxies are evolving on timescales that are significantly 
longer than this, our estimate should agree with those that 
use shorter-timescale indicators such as emission lines or IR 
photometry. In the next subsection, we will test this asser- 
tion using a catalogue of galaxy spectra generated using a 
cosmological simulation. 

In Figures 6 and 7, we show the probability distribu- 
tions of the parameter SSFR for each mass bin of the SDSS 
and DEEP2 galaxy samples. The solid lines show SSFRs 
estimated using BC03 models without dust, while dashed- 
lines show the results from models with dust. As can be seen, 
the results agree with each other to within 0.15dex, because 
dust has a very weak effect on the Balmer absorption line 
series. From now on, we will quote results based on models 
without dust. 

In Figure 8, we plot the median values of the proba- 
bility distribution functions (PDF) of SSFR as a function 
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Figure 8. SSFRs as a function of stellar mass. Green circles show 
the SDSS results. Blue triangles show the fibre SSFRs estimated 
by B04. Black diamonds show aperture-corrected SDSS SSFRs. 
The blue triangles and black diamonds have been slightly dis- 
placed in the x-axis (by minus and plus 0.05dex, respectively) to 
allow the error bars to be distinguished. Red asterisks: SSFRs for 
the four DEEP2 mass bins. 



of stellar mass for the galaxies in the SDSS (green circles) 
and DEEP2 (red asterisks) samples. The error bars on each 
point indicate the 68 percentile range in the PDF around the 
median. As can be seen the widths of the PDFs are quite 
narrow, indicating that the formal error on our estimates of 
SSFR are small. 



5.1 Aperture effects 

When we compare the low redshift and high redshift mea- 
surements we must account for aperture effects. 

The SDSS spectra are obtained through 3" diameter 
fibres. At the median redshift of our sample, only the cen- 
tral 25% of the galaxy light has been included, so aper- 
ture effects will affect the SSFR estimates, particularly for 
early-type spiral galaxies where the light from the central 
few kpc will be dominated by the bulge component of the 
galaxy and light from the disk will be underrepresented. The 
DEEP2 spectra are taken through 1" slits; the typical galaxy 
in our high- z sample has a diameter of ^ 2". Aperture ef- 
fects will thus be less of an issue here, unless SFR gradients 
are significantly stronger at high z, which seems unlikely 



l|Dominguez-Palmero fc Balcellsll2008h 



When we compare the SDSS and DEEP 2 SSFR mea- 
surements in subsequent sections, we must use aperture- 
corrected values for the SDSS galaxies. We do this by ap- 
plying the corrections estimated by B04. They use the g — r 
and r — i colors measured outside the fibre to estimate the 
amount of sta r formation in the o uter parts of the galaxies. 
We note that ISalim et all (|2007h showed that this correc- 
tion is an overestimate for certain classes of galaxies, such 
as AGN and galaxies without detectable Ha emission. Com- 
paring the average aperture-corrected SSFR of B04 in each 
mass bin with the results from Salim et al., we find that 
except for the most massive bin (log M*/M® = 11.5 — 12), 
the discrepancy is negligible. We have thus decided to ap- 
ply the B04 correction to our 6 SDSS mass bins and the 
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results are shown in Figure 8 as black diamonds (they have 
been displaced in the x-axis by plus 0.05dex to allow the 
error bars to be distinguished). The aperture correction in- 
creases the value of SSFR by ~ 0.15dex at M„ < 1O 1O M 
and 0.2 ~ 0.3dex at M» > 10 10 M Q We have not applied any 
similar correction to the high-z data. We will discuss this 
issue in more detail in §6.2. 

Comparing our results for SDSS and DEEP2 presented 
in Figure 8 leads to two main conclusions: 

(i) At stellar masses above I0 10 Mq, in both the low 
redshift and the high redshift samples, the SSFRs de- 
crease monotonically and relatively steeply with stellar mass 
(SSFRocM- (0 ' 85 ~ a9) ) . 

(ii) The average SSFR decreases by a factor of ~ 3 — 4 
from z ~ 1 to z ~ 0, and this factor appears to be indepen- 
dent of galaxy mass. 

5.2 Tests of the method 

We have claimed that our method of estimating the average 
SSFR of a stacked sample of galaxies using .ffbaimer should 
agree with other, more traditional methods that make use 
of indicators that are sensitive to star formation over much 
shorter timescales. 

In this section, we test this claim in two different ways. 
Our first test is a direct one. B04 have estimated SSFRs for 
the SDSS galaxies in our sample using emission line fluxes 
(the greatest weight is carried by Ha). For AGN and galaxies 
without detectable Ha emission, the SSFR is estimated from 
the 4000 A break strength measured within the fibre, using 
the relation between D4000 and SSFR calibrated using star 
forming galaxies. 

The blue triangles in Figure 8 show the average SSFRs 
from B04 for exactly the same sample of galaxies used to 
derive the green circles. We shift the blue triangles in the 
x-axis direction by minus 0.05dex to make the comparison 
clear. We find very good consistency between the two meth- 
ods, suggesting that .ffbaimer does provide a good measure 
of the average SSFR of a population of galaxies and that 
our methodology is robust. We note that the SSFRs de- 
rived by B04 from nebular emission lines agree well with 
those derived from multi-wavelength SED fitting including 
UV fluxes by Salim et al. (2007). This provides yet another 
confirmation of the robustness of the method. 

We cannot test our estimates of SSFR for the DEEP2 
galaxies in the same way, because Ha has shifted out of the 
relevant spectral range. Instead, we test our method using 
a DEEP2 mock catalogue of galaxy spectra that has been 
gener ated from the Millennium Run simulation (|Springell 
120051 ). We first construct a 2 x 2 sq. deg. moc k catalogue us- 
ing the MoMaF software l|Blaizot et al J 12003) and the semi- 
analyt ic prescription for galaxies from |Pe Lucia fc Blaizotl 
(2007) 0- An impo rtant ingredient of t his model, inherited 
from the work of ICroton et al.1 (|2006h . is the AGN feed- 
back implementation, which prevents cooling flows in mas- 
sive haloes, and hence quickly shuts down star formation in 
massive galaxies. In lower mass haloes, where cooling con- 
tinues unhampered by AGNs, the rate at which new cold gas 

3 The data for this model are publically available at 
http: / / www. mpa-garching. mpg.de / millennium 
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Figure 9. SSFRs as a function of -Hb a i mcr f° r model galaxies in 
the DEEP2 mock catalogue with redshifts in the range 0.75 ^ 
z $C 1. Blue lines mark the mean SSFRs of galaxies in the mock 
catalogue evaluated directly from the stellar masses and SFRs 
predicted by the semi-analytic model. Red lines mark the values 
of HbaLmei measured from the composite spectra. The green dots 
show the SSFRs derived from i?balmer- The last panel has so few 
objects that they are shown as black points rather than greyscale. 



reaches the galaxy and the efficiency of galactic winds reg- 
ulate the rate at which stars form. We select mock galaxies 
to match the DEEP2/POWIR magnitude selection, namely 
with 18.5 < Rab < 24.1 and A' AB < 22. We compute the 
stellar SED of each selected mock galaxy using its complex 
star formation history and metallicity evolution predicted 
by the semi-analytic model (SAM), combined with the BC03 
stellar population library. In the end, for each mock galaxy, 
we thus have an "observed" SED plus a set of properties (z, 
M», SFR) taken directly from the output of the SAM. Al- 
though this catalogue may not be a perfect representation of 
the real Universe, the co-moving star formation rate density 
does increase at higher redshifts at a rate similar to that 
observed (see Kitzbichler & White 2007 for a discussion). 
The catalogue should thus provide a good way to test our 
methodology using a sample of galaxies that is forming stars 
more actively than those in the local Universe. 

As seen in Figure 1 of W07, the strength of Balmer 
absorption lines decreases not only for old stellar popula- 
tions, but also for the very youngest stellar populations, 
which are dominated by O stars with weak absorption lines. 
If very young galaxies dominate the high redshift samples, 
then one might worry that the method might underestimate 
the specific star formation rates, simply because an O-star 
dominated stellar population would be mistaken for a much 
older one. We use the DEEP2 mock catalogue to show that 
at redshifts ~ 1, the average galaxy is not sufficiently young 
that this becomes a serious concern. 

We have selected galaxies with 0.75 ^ z ^ 1 from the 
DEEP2 mock catalogue and divided them into mass bins in 
exactly the same way as was done for the real data. The 
analysis of the model galaxies (including the stacking and 
filtering procedure) is also carried out in the same way as 
for the real data. In Figure 9, we plot SSFR as a function of 
#baimer for each individual model galaxy. Results are shown 
for four different mass bins. The blue lines show the average 
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SSFR of all the galaxies in the bin, defined as (SFR)/(M„). 
Red lines are the values of Hbalmei measured from the com- 
posite spectra. The green dots show the SSFRs derived from 
Maimer of composite spectra. As can be seen, the values ob- 
tained from the stacks are almost exactly the same as those 
obtained from the average instantaneous SSFR calculated 
from the model galaxies. 



6 DISCUSSION 

In this section, we translate SSFR into a dimensionless star 
formation activity parameter (see Dave 2008), defined as 



1 



SSFR t H (z) - lGyr' 



(15) 



Physically, this can roughly represent the fraction of the 
Hubble time (minus a Gyr) that a galaxy needs to have 
formed its stars at its current rate in order to produce its 
current stellar masfl A Gyr is subtracted in order to take 
account of the fact that dark matter halos massive enough 
to host galaxies with reasonably high star formation rates 
take about 1 Gyr to assemble in a ACDM Universe (iDavel 
2008). A value of a s f of 1 indicates that the galaxies could 
feasibly have formed all their stellar mass by forming stars 
continuously at the rate now observed. If a„f is greater than 
1, their past average SFR must have been greater than their 
current SFR for the stellar mass of the galaxy to have formed 
within a Hubble time. 

In Figure 10, we show our estimate of this characteristic 
timescale a s f as a function of stellar mass at both high (red 
asterisks) and low (black diamonds) redshifts. We remind 
the reader that the \ow-z, highest mass bin is uncertain due 
to aperture correction complications in the SDSS survey. 
From this plot, we find that a s f is a strongly increasing 
function of M, ; the timescale increases from values close to 
a Hubble time for galaxies with 1O 1O M0 to more than an 
order of magnitude larger than the Hubble time for galaxies 
more massive than 10 Mq. We find that, at fixed stellar 
mass, loga 3 / evolves only slightly with redshift, decreasing 
by around 0.2dex between z ~ and z ~ 1. 



6.1 Comparison with galaxy formation models 

The blue squares in Figure 10 (displaced in the x-axis by 
plus 0.05dex to make the comparison clear) show the char- 
acteristic timescale of star formation in the mock DEEP2 
universe created from the Millennium run SAM as described 
in Section 5.2. A key result of this paper is that the model is 
entirely consistent with observational results at z ~ 1, both 
in amplitude and slope. 

More generally, in theoretical models of galaxy forma- 
tion, at a fixed stellar mass, ot 3 f is predicted to r emain 
const ant out to redshifts greater than 2 (Figure 2 of lDavel 
2008). This latter result is robust to methodology, both semi- 
analytic and smooth-particle-hydrodynamic (SPH) simula- 
tions agree. Specifically, between a redshift of 1 and the 
models predict changes in loga a / of less than O.ldex. Our 

4 We note that a s f is not simply the fractional Hubble time the 
galaxy takes to form its stars, because stellar mass is returned to 
the ISM, a fact which is not accounted for in this simple model. 
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Figure 10. Star formation activity parameter a a f as a function 
of stellar mass. As in previous plots, black diamonds and red 
asterisks represent SDSS and DEEP2 results, respectively. Blue 
squares are the results derived from the DEEP2 mock catalogue, 
they have been displaced in the x-axis by plus 0.05dex to make 
the comparison clear. 



observational results indicate a slightly larger increase of 
0.2dex. Further studies using larger surveys will be required 
to confirm this small amount of evolution. 

The amplitude of a s f and its variation with stellar mass 
depend on simulation methodology. The near unity a s f at 
both z ~ and z ~ 1 for galaxies with M* < 10 ia5 M© 
suggests that galaxy mass growth in this mass and redshift 
range may be dominated by smooth and steady cold mode 
accretion, as implemented in all current models of galaxy 
formation. The strong increase of a s f with M* that we ob- 
serve at both redshift 1 and in the local Universe (Figure 
10), is an expression of the phenomenon of 'downsizing': 
massive galaxies have apparently completed most of their 
star formation at higher redshifts than low mass systems. In 
many current models of galaxy formation, the explanation 
of this behaviour is "AGN feedback" . More massive galaxies 
are more likely to host massive black holes which have the 
capability of producing more energy. Of equal importance, 
these massive galaxies are hosted by larger halos, where the 
AGN energy can be well coupled to the material that would 
otherwise cool and fuel star formation in the galaxies. As 
a result, AGN feedback through heating of the interstellar 
and intergalactic gas is more efficient in massive galaxies. 

In summary, the observed amplitude and evolution of 
OL s j as presented in this paper, provide firm constraints for 
all galaxy formation models. 



6.2 Comparison with previous results 

As we have shown, our results appear to be in good general 
agreement with cosmological galaxy formation models. We 
now turn to a comparison with previous results from the lit- 
erature. As we shall show, there is a significant discrepancy 
of 0.3 — 0.4dex between our -HWmcr derived SSFRs, and 
those derived primarily from multiwavelength broad band 
photometry. 

The results from two other studies of the SSFRs of 
galaxies at z ~ 1 are compared with our results in Fig- 
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Figure 11. SSFR as a function of stellar mass. Black diamonds: 
aperture-corrected SDSS data. Red asterisks: SSFRs of DEEP2 
galaxies. Blue crosses and green squares: SSFRs from Z07 & 
AEGIS respectively. 



ure 11. The blue crosses represent SSFRs derived from the 
COMBO-17 sample with Spitzer 24/xm and GALEX data 
(Zheng et al. 2007, hereafter Z07). These autho rs use the 
meas ured UV+IR luminosities to derive SFRs (jBell et al.l 
l2005h . They take account of the contribution from galaxies 
that are not individually detected at 24//m by stacking their 
images. 

The green squares show results calculated by us for a 
sample of galaxies selected from the All- Wavelength Ex- 
tended Groth Strip International Survey (AEGIS). The 
galaxies are a subsample of those used to create the com- 
posite spectra in this paper. This sample is not exactly the 
same as presented in Noeske et al. (2007, hereafter N07) and 
iNoeske et aT] (|2007bh . as only star-forming galaxies were in- 
cluded in their analysis. Instead, we have averaged the SS- 
FRs over all the galaxies in each stellar mass bin (i.e. both 
star- forming and quiescent ones), using the same weighting 
factors that we used in our spectral stacking analysis. Aside 
from this difference in sample, the method used to derive 
SFRs for the AEGIS galaxies is the same as used by N07. 

The method used by N07 to derive SFRs is slightly dif- 
ferent to that employed by Z07, in that information from 
emission lines in the spectra of the galaxies is utilized if it 
is available. For galaxies with /24 Mm > 60/xJy and strong 
emission lines, the total SFR is derived from a combination 
of the IR measurements and from DEEP2 emission lines 
(Ha, H/3, or [O ll]A3727, depending on z) with no extinc- 
tion correction. SFRs are derived from extinction-corrected 
emission lines only for blue galaxies with strong emission 
lines and no detectable 24 micron emission. Red galaxies 
with weak emission lines, but no 24 micron detections, are 
considered star-forming and SFRs are derived from emis- 
sion lines, assuming the same extinction corrections as for 
normal star-forming galaxies. 

The dashed lines in Figure 11 are a linear fit to the data 
points for mass bins with log M,/Mq ^ 10. As can be seen 
from this plot, the slope of the relation between SSFR and 
mass that we derive is consistent with the results of AEGIS 
and Z07. However our method yields a normalisation that is 
0.3 - 0.4dex lower than AEGIS and Z07. 



This discrepancy is puzzling, but it does have a number 
of possible explanations: 

1) Systematic differences in the calibration of different 
SFR indicators. Our method is based on stellar absorption 
line indicators, while the N07 and Z07 results are based on a 
combination of UV, IR and emission lines. We have demon- 
strated that our results do agree with the SSFRs derived by 
B04 from emission lines at low redshift, which provides con- 
fidence that there is no significant discrepancy between our 
method and the standard calibration of extinction-corrected 
Ha to derive SFR. We are unable to carry out the same test 
at z ~ 1, because Ha is redshifted out of the spectral range 
covered by most galaxy surveys. Systematic differences in 
estimated stellar mass do occur when different population 
synthesis models are used due to differing mass-to-light ra- 
tios in the models. However, the same BC03 models have 
been used in the comparison to the AEGIS results. 

2) Obscured AGN. There has been no attempt to re- 
move obscured AGN from the two 2 = 1 galaxy samples 
with which we compare. Both the N07 and Z07 analyses 
make use of the 24 micron Spitzer passband as a star for- 
mation indicator. At z = 1, this corresponds to a rest-frame 
wavelength of 12 micron, which is very close in wavelength to 
where emissio n from a dusty tor us would become very signif- 
icant (see e.g. lDaddi et al.ll2007M . In addition, AGN emission 
could well be contaminating some of the optical emission 
lines used to estimate SFR. By contrast, the -Hbaimcr index 
originates from stellar atmospheres and is not expected to 
be contaminated from emission from an obscured AGN. 

On the other hand, the vast majority of DEEP2 24 mi- 
cron sources and line-emitting galaxie s have line ratios in - 
dicating star formation and not AGN (jWeiner et al.l [20071 ). 
If the AGN were highly obscured, they could contribute at 
24 micron and not show up in optical lines, but to make 
up a difference of 0.3dex one would have to assign half of 
the 24/^m emission at z=l to obscured AGN, which would 
appear quite extreme (B. Weiner, private communication). 

3) Evolution of IMF with redshift. The SFR indicators 
used by N07 and Z07 trace O and B stars, whereas our 
SFR indicator is the Balmer series and is mainly influenced 
by A stars. If the IMF changes wit h redshift such that more 
mass ive stars form at higher z (e.g. |Dave1l2008l : Ivan Dokkuml 
2008), then a discrepancy between the two methods may not 
be apparent in the analysis of the low redshift samples, but 
may become more pronounced at higher z. 

4) Aperture effects. As we have discussed, the SSFRs we 
estimate for the z = 1 galaxies may be biased somewhat low 
because the long-slit spectra preferentially sample the inner 
bulge of the galaxy. However, 95% of the DEEP2 galaxies (at 
all redshifts) with a line measure ment have r e ff< .95" as 
measured in the CFHT imaging ijWeiner et alJl2007t ). Thus 
the 1" slit covers a large fraction of the galaxy. Additionally 
the seeing mixes the light into the slit to a much greater 
degree in DEEP 2 than it does for the SDSS fibers. Thus the 
star formation gradients of z = 1 galaxies would have to be 
extremely strong for this to make a factor 2 — 3 difference to 
the SSFR estimated for the galaxy population as a whole. 
As can be seen from Figure 8, the correction for aperture 
effects in the low z sample is only a factor of 2 on average, 
even in the case where the fibre only samples 25% of the 
total light. 

In this paper, we are not able make a definitive con- 
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elusion with regard to the possibilities listed above. It is 
clear that there are many inherent uncertainties in estimat- 
ing SFRs in galaxies and that more work is needed before 
the factors of 2 — 3 offsets that we see between the different 
methods and the models can be understood in detail. 



7 SUMMARY 

In this paper, we developed a new method to measure the 
average SSFR of a population of galaxies using the Balmer 
absorption line series located in the rest-frame wavelength 
range 3750 — 4I50A. Our method is free of complications 
due to dust extinction and AGN contamination, and pro- 
vides a consistent way to measure SSFRs at both high and 
low redshifts. The robustness of the method has been tested 
using SDSS data and a DEEP2 mock catalogue drawn from 
cosmological simulations. 

We apply this method to the DEEP2 galaxy sample. 
We high-pass filter the spectra because the DEEP2 spectra 
are not flux calibrated and then stack together the spectra 
of galaxies with similar stellar masses. SSFRs are estimated 
from the Balmer absorption line series by comparing to a 
library of model spectra generated from the BC03 stellar 
population code. 

Our results show that: 

• The average SSFR decreases monotonically with stel- 
lar mass at both z ~ and z ~ 1. At both redshifts the 
decrease is almost an order of magnitude between 1O 1O M0 
and 10 115 M©. 

• For galaxies of fixed stellar mass, the average SSFR has 
decreased by a factor of 3 — 4 from z ~ 1 to z ~ 0. 

• The amplitude of the decrease in SSFR with z is inde- 
pendent of stellar mass in the stellar mass range observable 
in the high redshift sample, i.e. 10 10 - 10 U,5 M Q . 

• The average SSFR of galaxies at z ~ I is consis- 
tent with the prediction s of the semi-analytic model of 
|Pe Lucia fc BlaizotJ (|2007T l in all stellar mass bins observed. 

We define a star formation activity parameter a s f which 
is the average time (in units of the Hubble time at redshift 
z) needed by a galaxy of a given mass to form its stars at 
its current SFR. We find that galaxies of mass < 1O 1O ' 5 M0 
are consistent with a s f ~ 1. This is in good agreement with 
models in which star formation is regulated by the infall 
of cold gas. At higher masses a s ; increases sharply, which, 
in the currently favoured galaxy formation scenario, can be 
understood by AGN feedback reducing the rate at which gas 
cools onto more massive galaxies. 
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